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Abstract We present models of temperature distribution in 
the crust of a neutron star in the presence of a strong toroidal 
component superposed to the poloidal component of the mag- 
netic field. The presence of such a toroidal field hinders heat 
flow toward the surface in a large part of the crust. As a 
result, the neutron star surface presents two warm regions 
surrounded by extended cold regions and has a thermal lu- 
minosity much lower than in the case the magnetic field is 
purely poloidal. We apply these models to calculate the ther- 
mal evolution of such neutron stars and show that the low- 
ered photon luminosity naturally extends their life-time as 
detectable thermal X-ray sources. 
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1 Introduction 

It is generally considered that, within a few decades after its 
birth in a core collapse supernova, a neutron star reaches a 
state of "isothermality" characterized by a uniform high in- 
terior temperature. The stellar surface is much colder and 
protected from the hot interior by a thin (~ 100 meters) 
layer, the envelope, which acts as a h eat blanket. As was 
argued by iGreenstein & Hartkel ([1983), in the presence of 
a sufficiently strong magnetic field, > 10 10 G, the surface 
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temperature of the neutron star will not be uniform as is 
expected in the unmagnetized case since the magnetic field 
severely limits the ability of electrons to transport heat in 
directions perpendicular to itself. As a result, the regions 
around the magnetic poles, where the magnetic field is al- 
most radial, are expected to be significantly warmer than 
the regions around the magnetic equator, where the field is 
almost tangent to the surface. Since then much work has 
been dedicated to study the effects of the magnetic field on 
the properties of the neutr on s tar envelope and cru st (see 
IPotekhin & Yakovlevll200ll and lPotekhin et alJl2003l for re- 
cent works). In the presence of a sufficiently strong mag- 
netic field, > 10 12 — 10 13 G, the anisotropy of heat trans- 
port actually extends to much higher densities and can even 
be present within the whol e crust. Recently, we have shown 
(Gepp ert et ai1l2004 l2006h that, in cases where the field ge- 
ometry in the crust is such that the meridional component of 
the field dominates over its radial component in a large part 
of the crust, the non-uniformity of the temperature, previ- 
ously considered to be restricted to the envelope, may actu- 
ally extend to the whole crust. The largest effect was found 
when a strong toroidal component was included in the crust, 
superposed to the poloidal component. This result, that the 
geometry of the magnetic field in the interior of the neutron 
star leaves an observable imprint at the surface, potentially 
allows us to study the internal structure of the magnetic field 
through modelling of the spectra and pulse profiles of ther- 
mally emitting neutron stars. 

There exist growing observational evidence that the an- 
isotropy of heat transport in the envelope alone, assuming an 
otherwise isothermal crust, can not explain the surface tem- 
perature distributions of some observed neutron stars. In the 
case of seve ral of the "Magnificent Seven" (see, e.g.,[Haberl 
|2004 l2007u optical broad band photometric detections can 
be interpreted as being due to the Rayleigh-Jeans tail of a 
blackbody. However, these optical data are well above the 
Rayleigh-Jeans tail of the blackbody detected in the X-ray 
("optical excess") and indicate the presence of an extended 
cold component of much larger area than the warm compo- 
nent observed in X-ray, the latter having an emitting radius 
(^3 — 5 km) much smaller than the usually assumed radius 
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of a neutron star (~ 10 - 15 kml. lSchwope et al.l d2005h tried 
to fit the lightcurve of RBS 1223 and concluded that only a 
surface temperature profile with relatively small, about 4- 
5 km s across, hot pola r regi ons may expla i n the observa- 
tions. IPons et all d2002h and iTriimper et all d2004h had ar- 
rived qualitatively at the same conclusion when they fitted 
the combined X-ray and optical spectrum of RX J1856.5- 
3754. In both cases, the smallness of the hot region is much 
below what can be reached by considering anisotropic heat 
transport limited to only a thin envelope. 

Little is known about the magnetic field structure in neu- 
tron stars which is very likely determined by processes dur- 
ing the proto-neutron star phase and/or in a relativel y short 



The source term H includes all possible "heating mecha- 
nisms" which will be neglected in the present work; how- 
ever the decay of a strong magnetic field ca n result in signif- 
icant heating dKaminker et al J 120061 120071) . Only the basic 
points will be presented here and w e ref er the reader to th e 
reviews 



Yakoylev & Pethickl (|2004|) and iPage et all d2006h . 
or lPage et al.1 ( 2004 ). for more details. 



The specific heat C v receives its dominant contribution 
from the baryons in the core of the star and about 10% con- 
tribution from the leptons. The crustal lattice, free neutrons 
in the inner crust and electrons also provide a small contri- 
bution. When nucleons become superfluid (neutrons) or su- 
perconductin g (protons) their contribution to C v is severely 
perio d after that epoch. A proto-neutron star dynamo (ThompsorediiBtuhaad may even practically vanish. 
19931) is unlikely to generate purely poloidal fields while dif- The neutrino luminosity L v is usually dominated by the 
ferential rotation will easily wrap any poloidal field and gen- core. We consider the slow emission processes of the modi- 
erate strong toroida l components (Klu zniak & Ruder man 1998;fied Urea fa mily and the similar bremsstrahlung ones. Once 
Wheel er et al .I2002h . The magneto-rotational instability (BalbuavfelBamtering turns on, neutrino emission from the forma- 
1991) also most certainly acts in proto-neutron stars dAkivama flt oih. and breaking of Cooper pairs is also included. 
2003 ) and results in toro idal fields from differential rotation 
(Balb us &Hawleyl 1998b . Thus, it seems realistic to consider 



the effect of magnetic field configurations which consist of 
poloidal and toroidal components. 

Besides their possible relevance for modelling the ob- 
served thermal radiation of the "Magnificent Seven", and 
probably other strong field isolated neutron stars, toroidal 
magnetic fields may have a strong effect on the thermal evo- 
lution of the stars. The highly non-uniform surface temper- 
atures they can induce result in reduced thermal luminosi- 
ties and hence reduced energy losses during the late photon 
cooling era. It is our purpose in this work to explore this is- 
sue and we will present preliminary models of cooling of 
neutron stars with strong toroidal fields. 

The structure of this paper is the following. In § [2] we 
briefly review the basic ingredients and concepts involved in 
modelling the cooling of a neutron star and in §[3] we present 
in some detail the simple mathematical formalism which de- 
scribes dipolar fields, both poloidal and toroidal. The next 
section, § [4] presents our results on the effect of strong mag- 
netic fields on the neutron star crust temperature distribution, 
and these results are applied to model the cooling of the star 
in § El Finally, we discuss our results in § [6] and offer some 
tentative conclusions and prospects for future work. 



2 Basic mechanisms of neutron star cooling 

The essentials of neutron star cooling are expressed in the 
energy balance equation. In its Newtonian formulation this 
equation reads 



2.1 Neutron star envelopes (without magnetic fields) 

Our main interest here is the photon luminosity Ly. In the 
absence of a magnetic field one expects the surface to have 
a uniform temperature T e , the effective temperature, and can 
express Ly as 

L y = AnR 2 OsbT* or L yoo = AnRl a SB T^ (2) 

where R is the star's radius and Osb the Stefan-Boltzmann 
constant (the quantities "at infinity" are defined as Ly„ = 
e 20 L r , = e- R, and T eoa = e T e , where e = (1 -2GM/Rc 2 ) 
is the redshift factor). In order to integrate Eq. [T|one needs 
a relationship between T and T e . The assumption of uni- 
form interior temperature T is reasonable for most of the 
interior, given the huge thermal conductivity of degenerate 
matter and once the star is old enough to have relaxed from 
the initially complicated temperature structure produced at 
its birth, but is certainly not possible in the upper layers of 
the star where density is low enough for the electrons not 
to be fully degenerate. One traditionally separates out these 
upper layers from the cooling calculation and treats them 
as an envelope. A typical cut density is 10 10 g cm -3 and 
the resulting envelope, with a depth of the order of 100 me- 
ters, can be studied separa t ely in a pla ne parallel approxima- 
tion. Gundmundsson et ail (Il982l) (and lGundmundsson et al.l 
Il983l) presented detailed models of neutron star envelopes 
and showed that T e is related to the temperature at the bot- 
tom of the envelope, Tj, at density Pb 
the simple relation 
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10 10 g cm 3 , through 



dt 



dT 
dt 



-L v — Ly + H . 



(1) 



71 = 0.87 x 



io 6 k^: 



-0.55 



(3) 



where is the thermal energy content of the star, C v its 
total specific heat and T its internal temperature, which we 
assume to be uniform here. The energy sinks are the neu- 
trino luminosity L v and the surface photon luminosity Ly. 



where g s u is the surface gravity in units of 10 cm s . 
A relation such as Eq. [3] is usually called a "Tj, — T s re- 
lationship". The models leading to Eq. [3] assumed the en- 
velope is formed of iron-like nuclei, and it was sho wn by 
IChabrier et all dl997h (see also iPotekhin et alJll997h that if 
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Fig. 1 Magnetization parameter a> B z vs. density at six different tem- 
peratures (as labeled on the curves) assuming a uniform magnetic field 
of strength B = 10 12 G. I ts value for d ifferent field strengths scales 
linearly in B. (Figure from Geppert et al. 20041.) 



light elements, such as H, He, C, O, are present deep enough 
in the envelope, the increase in thermal conductivity (which 
is roughly proportional to Z _1 in liquid matter, Z being the 
element's charge) results in much higher luminosities, by up 
to one order of magnitude. Since the temperature gradient 
penetrates deeper in hotter stars, larger amounts of light ele- 
ments are necessary to alter heat transport at high Tj, than at 
lower 7],. 



2.2 Heat transport with magnetic fields 

In the absence of a magnetic field the thermal conductivity 
K is conveniently written as 



1 -2 

Kb = -C V V Z X 



7Z 2 kftTn e 

3m* 



(4) 



where (Ob = eB/m*c is the electron cyclotron frequency. The 
condition (ObT ^> 1, which implies strong anisotropy, is eas- 
ily realized in a neutron star envelo pe, and also possibly in 
the whole crust (Geppert et alj2004l) . Values of the magneti- 
zation parameter (ObT are plotted in Fig.Q] Notice that T re- 
ceives contribution from both electron-phonon and electron- 
impurity scattering, with frequencies V e - P h and V e _i mp resp., 
and is given by T = (v e - p h + V e _i mp ) _1 . Electron-impurity 
scattering dominates at low temperatures, and/or high densi- 
ties, and is T-independent while V e - P h goes roughly as T 2 : 
combination of these two different T -dependencies is the 
reason for the complex behavior of (OgZ seen in Fig.Q] 

In case the field is strong enough to be quantizing, the 
expressions [6] have to be modified (in particular T also be- 
comes anisotropic) but t he essential re sult that K± <C Kh, 
when ObT 3> 1, remains dPotekhinlll999h . 



2.3 Magnetized neutron star envelopes 

In considering magnetic field effects on heat transport in a 
neutron star the simplest case to handle is the envelope: since 
its thickness, ~ 100 m, is much smaller than the length scale 
over which the field is expect to vary significantly, ~ km, 
one can consider heat transport on a given small patch on the 
surface to be independent of the rest of the surface. More- 
over, the field strength B and its angle with respect to the ra- 
dial direction, 0g, can be considered as uniform in the patch 
and the heat flux considered as essentially radial. The ther- 
mal conductivity in a direction making an angle 0b with the 
field, considering the jc tensor of Eq.|5] is then given by 



k{&b) = cos 2 @ B x K|i + sin 2 Q B x k 



(7) 



where v is the mean velocity of the heat carriers, c v their spe- 
cific heat per unit volume, and T their collisional time; the 
second expression is particularized to relativis tic electrons, 
the do minant heat carriers in neutron star crusts (iYakovlev & U 
1980). In the presence of a magnetic field, due to the classi- 
cal Larmor rotation of electrons, heat flow may be anisotropic 
and K becomes a tensor 



With this form of k(&b) and a radial flux, heat transport in 
the envelope at this surface patch is a one dimensional prob- 
lem. The solution is a "7}, — T e " relationship which depends 
on B and 0g, and the obtained effective temperature T e is a 
local one which we will write as T s (6,<j>) in the sense that 
the flux emergin g from this patch is OftT*. G iven the form 
of k(&b), Eq. 171 iGreenstein & Hartkd d 19831) proposed the 
simple interpolation formula 



K = 




</>) 4 = 7;(r b ;B,0 B ) 4 « 

cos 2 (0 B ) x T s {T h ;B,0 B = O) 4 + 
sin 2 (0 B ) x r s (7 b ;S,0 B = 9O°) 4 



(8) 



(5) 



(assuming the field B oriented along the z-axis) whose com- 
ponents have the form 



K|| = Kb 



Kb 



1 + (<%T) 2 
Kb (ObT 



(6) 



1 



{(0 B T) 2 



for arbitrary angle 0b in terms of the two case s of radial 

B =| 0) an d tangential (0r = 90°) field. Recently. IPotekhin & Yakovlevl 
200 ll) and IPotekhin et ail d2003l) have presented detailed 

calculations and fitted their results by an expression similar 

toEq.[U 

For a given geometry of the magnetic field (not necessar- 
ily dipolar) and assuming that the temperature at the bottom 
of the envelope, 7],, is not affected by the magnetic field and 
hence is uniform arround the star, one can generate the ex- 
pected surface temperature distribution 7^ ( 9 , ) by piecing 
together envelope models through applying relationships of 
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Fig. 2 Two examples of surface temperature distribution obtained from Eq.[8] The maps cover the whole neutron star surface in an area preserving 
projection. For better viewing, the magnetic field symmetry axis is located in the equatorial plane, oriented from (j) = 90° to = 27 0°. The left 
panel has a dipolar field only while the right panel also contains a quadrupolar field superposed to the same dipole. (Figure from Pag e et all 
2006.) 



the form of Eq. [8l at eac h point of the surface (Page Hl995l : 
iPage & Sarmientol [l996). Two examples of such expected 
surface temperature distributions are shown in Fig. |2j The 
photon luminosity is then given by 



L y = 4nR z 



^o B T s (6,$) 4 ^4xR 2 
4% 



(9) 



where T e is again defined as in Eq.|2] 

By this mehod one can easily generate surface tempera- 
ture distributions corresponding to the geometry of the mag- 
netic field at the surface (assuming that the same field geom- 
etry is maintained throughout the underlying thin envelope). 
However, the results of Fig. Q] show that the anisotropy in 
heat transport is likely to extend much deeper into the crust 
than just the envelope. To handle such cases one is hence 
forced to define the field geometry in the whole crust, which 
is what we do in the next section. 



3 The internal magnetic field 

We will only consider here axisymmetric configurations (mo- 
re general ones will hopefully be considered in the future 
once a 3D transport code is available). In this case it is con- 
venient to separate the magnetic field in two components 



B = B pol +B t 



(10) 



the poloidal and toroidal components, respectively, where 
BP°i on iy has e r and eg components and B tor only an e^ 
component] (the e r e being units vectors in the spherical 
coordinate directions, with the 0=0 axis coinciding with 
the field symmetry axis). Thus, the magnetic field lines of 
B tor are simply circles centered on the symmetry axis but 
the field lines of B po1 are more complicated. Let us expand 
B r in Legendre polynomials as B r = £F/(r) P/(cos0). We 
will only consider the first (/ = 1) term, the dipole, and write 
F\ as S/r 2 so that 



Br 



S(r) 



cos 8 



(11) 



1 Such decomposi tion is also p ossible for non-axisymmetric fields 
but is more involved (Radlei 2000). 




Fig. 3 The three components of the magnetic field considered in this 
work 



Then Maxwell's equation V • B = implies 

1 dS . „ 
B e = - — ^-sinfj . 
2r or 

We also need the three boundary conditions 

S(r = 0) = 



S(r = R) = B R 2 and 



dS 
dr 



r=R 



R 



(12) 



-B Q R (13) 



which ensure regularity at the star's center and smooth match- 
ing with an external vacuum dipolar field, of strength So 
at the magnetic poles, at the stellar surface. Besides these 
boundary conditions, the Stoke function S is totally arbitrary 
but a choice of it is equivalent to choosing the location of the 
currents sustaining the poloidal field since the latter, having 
only a j§ component, are given by 

c sin0 fd 2 S 2S\ 

" = 5?— U?-?j (14) 

Notice also a simple physical interpretation of the Stoke func- 
tion: the magnetic flux through the star's equatorial plane in 
a circle of radius r is simply given by nS(r), and the bound- 
ary condition gives a total flux 4> = nR 2 Bo, as it should be. 
In vacuum, S simply reduces to BqR 3 jr. 

The locations of the currents in the stellar interior is un- 
known but it is natural to separate them in two components, 
located in the core and in the crust and accordingly separate 
BP Ql as 

gpol =B core +B crust _ ( , 5) 
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Fig. 4 Continuous lines: normalized Stokes function s(x) = S/R 2 B^" sl 
and its derivative (scaled by a factor 10) used in this work, which gener- 
ate B,. (Eq.QTJi and Bg (Eg. 1121 resp. Discontinuous lines: the three dif- 
ferent normalized functions t(x) = T /Bg or we consider for the toroidal 
field, Bp (Eq.[T6t. which w e label as Tl, T2, an d T3. On the horizontal 
axis x = r/R. (Figure from Geppert et al. 200^.) 

The core protons are expected to become a supe rconduc- 
tor, w ith critical temperatures T c ~ 10 9 K (see, e.g. JPage et alJ 
2004), soon after the star's birth: the magnetic field is then 
confined into flux tubes and maintained by proton supercur- 
rents. Since, by definition, B core corresponds to currents lo- 
cated in the core, within the crust it is described by a vacuum 
dipolar poloidal field, and we will parametrize it by its sur- 
face strength at the magnetic pole #Q 0re . We do not need to 
specify the distribution of the supercurrents sustaining B eore 
and the geometry of this component in the core since the 
field, being confined to fluxoids which occupy only a very 
small volume, is not expected to alter heat transport in the 
core. 

For the B crust component, we need to specify its geome- 
try in the crust, i.e., the corresponding Stoke function S crust . 
The arbitrariness involved in such specification can be some- 
what relieved by considering models of the time evolution 
of this component: currents spontaneously migrate toward 
the highest density region of the crust, where the electrical 
resistivity is smallest, until they reach the crust-core bound- 
ary where their migration is stopped by the p roton super- 
conductor (see, e.g., Fig. 4 in lPage et al.ll2000h . We use an 
S CTUSt function resulting from such evolutionary calculations 
and scale it to vary the overall strength of B crust , which we 
parametrize by Bp™ 1 *' defined as its strength at the magnetic 
pole. 

For the toroidal component B tor , we also expand it in 
Legendre polinomials P; (sin d ) and keep only the I = 1 dipo- 
lar term 

B tor = r(r)sin0. (16) 

We do not have to consider the part of B tor confined to the 
core, because of proton superconductivity, and specify only 
the part confined to the crust. The only restrictions are the 
boundary conditions 

T(r = R COYe )=0 and T{r = R)=0. (17) 




Fig. 5 Temperature distribution in a strongly magnetized neutron star 
crust (whose thickness has been stretched by a factor five for easier 
reading). The choosen field scale parameters are Bq 01 " 6 = 7.5x10 13 G, 
B£ n,st = 2.5 x 10 13 G, Bg"" = 3 x 10 15 G, and the toroidal component's 
generating functions T is the model "Tl" of Fig. |4j- The color code 
maps the relative temperature, i.e., T(r, 9)/T core , with a core tempera- 
ture T core = 6 x 10 7 K. White lines show field lines of BP o1 , the field 
lines of B tor being perpendicular to the plane of the figure. The heat 
blanketing effect of the toroidal component is clearly visible. (From 
IGeppert et all20M ) 

We are not aware of evolutionary calculations of toroidal 
field in neutron star crust and hence are left with a guess 
about t he possible shape an d size of the T function: fol- 
lowing Geppert et al] d2006l) we see it as a free parameter 
and consider several choices. We parametrize the strength of 
B tor , for each choice of T, by Bq ot defined as the maximum 
val ue of | |B tor | | in the crust. 

iPerez-Azorin et al.l ((2006) have presented similar mod- 
els of field structure and the resulting temperature distribu- 
tions in the crust and chose force-free field configurations 
to specify the field geometry. It is reasonable to assume the 
field could evolve into a force-free configuration during the 
early neutron star life, but there is no physical reason why 
the later evolution, driven by Ohmic decay and Hall inter- 
actions, will conserve the force-free condition and this mo- 
tivates our considering the function T as rather arbitrary. 
However, our choice of S cmst , inspite of being based on evo- 
lutionary calculations of purely poloidal fields, may not be 
too realistic because of the coupling between poloidal and 
toroidal components from the Hall term. We show in Fig. [3] 
a sketch of the field geometry we consider and in Fig.|4]our 
choices of the functions S and T. 



4 Magnetic field effects in the crust 

Given the mag netic field geometri es described in the previ- 
ous sections, IGeppert et alJ (|2006) calculated the resulting 
temperature distributions in the neutron star crust, in a sta- 
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tionary state. In theses calculations, heat transport is solved 
in the crust for densities from p = p core = 1 .6 x 10 14 g cm -3 
down to p = p b = 10 10 g cm -3 . In the core, at p > p cor& , 
the temperature r core is assumed to be constant and uniform 
while at the outer zone, at p = Pb, the models are matched 
with magnetized envelopes (treated as the outer boundary 
condition). Since the 2D transport code we use does not yet 
include correct specific heat it does not have the capabil- 
ity to perform realistc time dependent calculations and only 
stationary results can be obtained, i.e., the thermal evolution 
is followed until the temperature profile would not change 
anymore with time (this stationary limit is independent of 
the specific heat). Moreover, neutrino energy losses are also 
not included (neutrino emission, for high temperature, and 
also strong magnetic fields, is not negligible and should be 
includ ed for accurate calculations, see, e.g., Potekhin et al. 
l2007h . Fig. |5] shows an example of the resulting tempera- 
ture distribution in the crust and Fig. [6] shows the resulting 
surface temperature distribution: this latter figure should be 
compared with the surface temperature distribution shown 
in Fig. |2] where isothermal crusts were considered and mag- 
netic field affected heat transport only within the thin enve- 
lope. 

For a fixed magnetic field configuration and several val- 
ues of r coie , we calculate the crustal temperature distribution 
and the resulting surface temperature T s (6), which is <j> in- 
dependent because of axial symmetry, and obtain T e from 
Eq. [9] A set of results is displayed in Fig. [7] we show mod- 
els which, within our selection of field configurations, maxi- 
mize the magnetic field effects, with 75% of the poloidal flux 
coming from crustal currents and the toroidal component lo- 
cated in the middle of the crust. With a purely poloidal field, 
i.e., Bq ot = 0, one sees little effect and the crust stays close 
to isothermality, except at high temperature, T core = 10 9 K, 
where the envelope's temperature gradient extends deeper 
into the crust than the p = 10 10 g cm -3 arbitrary cut. Sig- 
nificant effects appear when B 1 ™ 3> 10 14 G and when Bq y > 
10 15 G the crust is highly non-isothermal even at the high- 
est temperature r core = 10 9 since then COqT 1 in the whole 
crust. There is an intriguing change in the shape of the T- 
surface when going from r core > 10 8 K to r core < 10 8 K, 
which is most evident at the highest value Bq ot = 3 x 10 15 
G: it is most probably due to the shift in the dominant scat- 
tering process from electron-phonon, at high T where T is 
T -dependent, to electron-impurity, at low T where T is T- 
independent, as mentioned in Section |4] (We hope to study 
this effect in more detail in a future work.). 



5 Cooling of strongly magnetized neutron stars 

In order to perform (time dependent) cooling calculations 
we will here use the models described in the previous sec- 
tion to produce a set of o uter boundary co nditions for our ID 
cooling code (see, e.g., Page et alj|2004l) . We fix the outer 
boundary at a density Pb = 10 14 g cm -3 , instead of 10 10 g 
cm -3 when envelope models are used as boundary condi- 




Fig. 6 Surface temperature distribution corresponding to the crust 
temperature shown in Fig. [5] For better viewing the magnetic sym- 
metry axis is in th e equator (9 = 90°) pointing at <j) = 90°. (From 
Geppert et al. 2006.) 



tions. As is obvious from Fig. [6] compared to Fig. [2] the 
photon luminosity, and hence T e , is much lower in presence 
of a strong toroidal field than when the crust is considered 
as isothermal and Fig. [8] shows the resulting Tj, — T e rela- 
tionships obtained for nine different magnetic field config- 
urations. For the most extreme case, T e can be reduced by 
a factor of 2.5, and hence Ly by a factor 40, compared to 
the isothermal crust case. Notice that we chose Pb = 10 14 g 
cm -3 inspite of our crustal models starting at p core = 1.6 x 
10 14 g cm -3 because we found a very small temperature gra- 
dient in the range 1.0 — 1.6 x 10 14 g cm~ 3 and hence prefer 
to treat this density range with the cooling code. The ID 
cooling code solves the energy balance and heat transport 
equations in their full general relativistic forms. 

The cooling models we wil l consider here are within the 
minimal cooling paradigm of iPage et al.l d2004l) . In short, 
this means neutrino emission in the stellar core is from the 
modified Urea and the similar nucleon bremsstrahlung pro- 
cesses, neutron and proton pairing is taken into account with 
the resulting neutrino emission from the Cooper pair break- 
ing and formation process and the alteration of the specific 
heat. We use pair ing critical tempeature s for proton 'Spfrom 
lTakatsukaldl973h . for neutron l S from lSchwenk et alJd2003h . 
and for neutron 3 / > 2 the model "a" from IPage et al.l d2004l) . 
The star model is a 1.4 neutro n star built with the equa- 
tion of state o f lAkmal et all d!998h and hence, as part of the 
minimal cooling paradigm, charged meson condensates^, hy- 
perons, and deconfined quark matter are not present in the 
star. 

First we assess the effect of the approximation of trun- 
cating the star at Pb = 10 14 g cm" 3 , and treating most of the 
crust through our stationary solutions as part of the outer 
boundary condition, instead of using the traditional outer 
boundary at 10 10 g cm -3 . A comparison of cooling trajecto- 
ries with these two different boundary densities, and in the 
absence of a toroidal field component, is shown in Fig. [9] 
The major difference appears at early times during the crust 
relaxation era: a "full" (i.e., using Pb = 10 10 g cm -3 ) model 
shows strong r adial temperature ra dients in the curst at this 
stage (see, e.g., Gnedin e t al .11200 lh while the truncated (i.e., 
using Pb = 10 14 g cm~ 3 ) model assumes an isothermal crust. 



2 The APR equation of state shows the presence of a 7t° condensate, 
but it has no noticeable effect on the cooling. 
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Fig. 7 3D plots of the temperature distribution in a strongly magnetized neutron star crust. The horizontal axes show radial coordinate x = r/R 
and polar angle 6 while the vertical scale is T/T COK . Four values of T core are considered and the choosen field scale parameters for the poloidal 
field are Bg° le = 2.5 x 10 13 G and Bg rusl = 7.5 x 10 13 G in all cases while the toroidal component scale B'™ is varied, the overall shape of the T 
function, Eq.[l6] being the model "T2" of Fig.0 



Another way of seing it is that our 2D crustal models only 
considered stationary states and are hence only applicable 
when the cooling time scale is much longer than the thermal 
relaxation time of the crust, a condition which is certainly 
not fulfilled during this early cooling phase. Later, during the 
neutrino cooling era, the truncated model is slightly warmer 
than the "full" model because neutrino emission from its 
truncated crust is missing. During the photon cooling era the 
difference between the two models becomes larger, now due 
to the smaller specific heat of the truncated model which 
consequently cools faster. Overall, moving the outer bound- 
ary from 10 10 g cm -3 to 10 14 g cm -3 only has a small, and 



quite negligible, effect except at early ages which we will 
hence not show in our next results. 

Having now the confidence that stripping the star of most 
of its crust, and properly including the stripped part into the 
outer boundary condition, introduces an acceptably small 
error, we can proceed with models having strong toroidal 
crustal fields. We show in Fig. [10] our results for the nine 
field configurations of the 7}, — T e relationships presented in 
the Fig. [8] During the neutrino cooling era, the star's core 
evolution is driven by L v and the surface temperature, and 
Ly, simply follows the evolution of the core. So during this 
phase all nine models have exactly the same evolution but 
they look very different at the surface: models with a higher 
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Fig. 8 Possible Tj, — T e relationships for a neutron star with a strong 
toroidal magnetic field. The various sets of curves correspond to two 
different values of Bg or and two different locations of the toroidal field, 
labelled as "Tl" and "T2" as in Fig.[4] For each case we also consider 
two different splitting of the poloidal field: B c mst = 7.5 x 10 12 G and 
B= ore = 2.5 x 10 12 G, or B= rusl = 2.5 x 10 12 G and B™ le = 7.5 x 10 12 
G, the cases with the larger Bq™ 6 resulting in higher T e . The continuous 
cuve shows the T\, — T e relationships for an isothermal crust with the 
same poloidal field Bp° l = 10 13 G and magnetic field effects included 
only in the envelope (Potekhin & Yakovlev 200T1). 



r e for a given 7b, as shown in Fig. have a higher photon 
luminosity. During the photon cooling era the results are in- 
verted since Ly drives the cooling and models with a higher 
T e for a given 7], result in a larger Ly and consequently they 
undergo faster cooling. 

We have preferred to plot the cooling curves as Ly vs t 
instead of T e vs t: given the highly non-uniform T s (0,(j)) the 
effective temperature T e looses any observational meaning. 
More detailed models t aking this into ac count will be pre- 
sented in a future work (Page et al. 20Q3). 




Log t [yrs] 

Fig. 9 Comparison of cooling trajectories, i.e., red-shifted surface pho- 
ton luminosity vs time, obtained with our truncated models using Pb = 
10 14 g cm~ 3 (continuous curve) vs. a "full" model having p b = 10 10 g 
cm~ 3 (dotted curve). Both models consider a poloidal dipolar field of 
strength B[j o1 = 10 13 G. The truncated model uses the 7b — T e relation- 
ship for an isothermal crust shown in Fig. [§] while the full model uses 
directly the Tb — T e relationship from Potekiiin & Yakovlev (2001). 
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6 Discussion and Conclusions 

We have extended our previous results about tempe rature une-sty 
distribution in a strongly magnetized neutron star crust dGeppert et al.l 
12004 120061) to a broader range of temperatures and applied 
them to study th e cooling of neutron s tars with strong toroidal 
magnetic fields. Ic-eppert et ail d2004l) had shown that with a 
purely poloidal field entirely confined to the crust strongly 
non-u niform tempe r ature distributions develop in the crust. 
Then lGeppert et al.1 (120061) showed that allowing part of this 
poloidal field to permeate the core significantly reduced the 
temperature non-uniformity but that the inclusion of a strong 
toroidal field component can result in crustal temperature 
gradients even stronger than in the first case (sim ilar results 
have been obtained by IPerez-Azorin et al.ll2006h . We have 
shown here that, due to the strength of the field, the magne- 
tization parameter (0 B z is very large in most of the crust even 
at core temperatures as high as 10 9 K and the large crustal 
temperature gradients are still present in such hot stars. Con- 
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Fig. 10 Cooling of neutron stars with strong toroidal magnetic fields: 
red-shifted surface photon luminosity vs time. Each curve corresponds 
to one of the 7b — T e relationships displayed in Fig. [8] with the same 
line-style. 



sidering these results we have been able to perform cool- 
ing calculations of such neutron stars with a strong toroidal 
crustal magnetic field, with some approximations which, as 
we have shown, introduce only very small errors in the re- 
sults. 

Our final results are displayed in Fig. [10] and compared 
with the estimated luminosity range of the "Magnificent Sev- 
en" ( IHaberlll2004l) . This comparison indicates most probable 
ages between 0.5 and 1.5 Myrs for these stars, depending 
on the magnetic field structure, but this range could be ex- 
tended down to 0.05 Myrs for the brightest ones and up to 3 
Myrs for the dimest ones if the most extreme field geometry 
is considered. Besides the uncertainty in the field structure, 
which in itself results in the range of predicted ages shown 
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in Fig. [lOl there are two other key ingredients in cooling 
models which can have similar effects: nucleon pairing in 
the core and the chemical composition of the envelope. The 
former significantly affects the specific heat during the pho- 
ton cooling era and different assumptions about the values 
of T c can introduce a factor of a few in the predicted cooling 
ages while the latter can significantly affect the photon lumi- 
nosity (see discussion in ij l2.ll) and introduce another uncer- 
tainty of a few (we refer the reader to the detailled presen- 
tation of lPage et al.ll2004h . Considering these three theoreti- 
cal uncertainties, magnetic field geometry, nucleon pairing, 
and envelope chemical composition, it is certainly possible 
to extend the theoretical cooling ages, for photon luminosi- 
ties between ~ 3 x 10 30 and ~ 3 x 10 31 , to cover the range 
from 10 5 to almost 10 7 years ! 

Can we reasonably expect to reduce this enormous the- 
oretic aj_a^ ur^rtainty? Precession of some neutron stars 
(lLinkll2007l : IPoTrs]l2 0Q7) may imply that n eutron pairi ng crit- 
ical temperatures in the core are very low dLinkl20 03), a pos- 
sibility which moreover h as a strong theoretical fundament 
dSchwenk & Frimanl2004h : this would imply a large specific 
heat and favor long cooling times (within the minimal cool- 
ing paradigm). Interpretation of the absorption lines detected 
in most of the "Magnificent Seven " (see, e.g.. lHaberll2007l 
and Ivan Kerkwijk & Kaplan! l2007h may provide crucial in- 
formation about the surface chemical composition. How- 
ever, atomic and molecular p hysics in stron g field still has 
many secrets to be unveiled (ITurbinen l2007) and, moreover, 
what is needed is the chemical composition of the deeper 
layers of the envelope, i.e., several tens of meters below the 
surface. Finally, the structure and evolution of the magnetic 
field with toroidal components, in the neutron star context as 
considered in this paper, is an almost uncharted territory and 
much progress can be expected. 
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